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ABSTRACT 

IceCube has reported the detection of a diffuse TeV-PeV neutrino emission, for 
which the flat spectrum radio quasars (FSRQs) have been proposed to be the can¬ 
didate sources. Here we assume that the neutrino flux from FSRQs is proportional 
to their gamma-ray ones, and obtain the gamma-ray/neutrino flux ratio by the diffuse 
gamma-ray flux from Fermi- LAT measurement of FSRQs and the diffuse neutrino flux 
detected by IceCube. We apply this ratio to individual FSRQs and hence predict their 
neutrino flux. We find that a large fraction of candidate FSRQs from the northern sky 
in the IceCube point source search has predicted neutrino flux above the IceCube up¬ 
per limit; and for the sample of stacking search for neutrinos by IceCube, the predicted 
stacked flux is even larger than the upper limit of stacked flux by orders of magnitude. 
Therefore the IceCube limit from stacking searches, combined with the Fermi- LAT 
observations, already rejects FSRQs as the main sources of IceCube-detected diffuse 
neutrinos: FSRQs can only account for < 10% (< 4%) of the IceCube-detected dif¬ 
fuse neutrino flux, according to the stacking searches from the whole (northern) sky. 
The derived neutrino/gamma-ray flux ratio is 1, also implying that the gamma-ray 
emission from FSRQs cannot be produced by the secondary leptons and photons from 
the pion production processes. The caveat in the assumptions is discussed. 

Subject headings: neutrino - gamma-ray: general - galaxies: active 


1. Introduction 


The detection of extra terrestrial diffuse TeV-PeV neutrino emission has been reported recently 


by IceCube collaboration (lAartsen et alJl2013c IlceCube Colla boration 


with an excess to the atmospheric background up to 5.7a (lAartsen et al. 


neutrino spectrum can be fitted with a flat form (Aartsen et al. 2014a) 


2013 


2014a). 


Aartsen et ah 2014al) 


The single flavor 


El$ Vi ic = (0.95 ± 0.3) x 10 _8 GeVcm _2 s _1 sr _1 


( 1 ) 






















- 2 - 


from tens TeV up to 2 PeV, assuming flavor ratio of 1:1:1. The lack of neutrino events above 2 
PeV may suggest a cutoff at the neutrino spectrum. The sky distribution of the neutrino events 


is consistent with being isotropic ( Aartsen et al. [2014al) and hence of extragalactic origin, but 


the low statistics so far cannot make any solid conclusion. There are many models proposed 
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Liu et al 

2013; 
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3 2006 
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and by 


cosmog enic neutrinos from u l trahigh energy CR interactions with cosmic microwave background 
(CMB) ( Kalashev et al. 2013 : Laha et al. 2013 : Roulet et al. 2013 ). By u sing the Fermi-LA T ob¬ 
servations and considering the gamma-ray and neutrino flux correlation, Wang et al. (2014) sug¬ 
gested that only starburst galaxies and AGN jets in the above mentioned scenarios may produce 
the neutrino flux detected by IceCube. We further investigate the AGN jet origin in this paper. 

AGNs have long been considered as u ltrahigh energy CR s o urces, and hen c e prod ucing neu¬ 
trinos through photo-pion processes. Ref Murase et al. (2014); Dermeretal. (2014) suggested 
that CRs accelerated in the inner jets of AGNs, especially flat spectrum radio quasars (FSRQs), 
interact with photons from the broad-line region (BLR), can produce neutrinos that may account 
for the IceCube detected PeV neutrino flux. The Fermi Large Area Telescope (LAT) already pro¬ 
vides us a deep survey of the sky from about 30 MeV to several hundred GeV. Based on the 2-year 
observation of Fermi-LAT, the FSRQs’ luminosity function (LF) and redshift distribution are well 
determined (1 A iello et al. 2012). Assuming a correlation between the gamma-ray and neutrino 
fluxes from FSRQs, we constrain the neutrino emission by FSRQs using gamma-ray flux from 
Fermi-LAT observations, in order to examine the FSRQ scenario for IceCube detected neutrinos. 
In Section 2 we describe the assumption and the method of the analysis, Section 3 shows the data 
and results, and conclusion and discussion are presented in Section 4. 


2. Assumption and method 


We assume that (1) the gamma-ray and neutrino fluxes from FSRQs are proportional to each 
other; and (2) FSRQs can account for the diffuse neutrino flux detected by IceCube. The following 
analyses are all carried out under these two assumptions. 


There are leptonic and hadronic models for blazar gamma-ray emission ( Bottcher et al.ll201~3 ). 
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In the leptonic model the FSRQ gamma-ray emission is produced by electrons accelerated in the 
inner jet upscattering background soft photons from self produced synchrotron radiation or the ac¬ 
cretion disk and BLR. On the other hand, in the hadronic model the FSRQ gamma-ray emission is 
mainly produced by CR proton synchrotron radiatioifl In both cases assumption (1) is reasonable 
if the energy ratio between the jet accelerated electrons and CRs is constant. Moreover, assumption 
( 1 ) should be more reasonable and valid in a statistic sense, although it may not tightly hold for 
individual FSRQs. 


The FSRQs are bright sources in the gamma-ray sky. The luminosity function (LF) and red- 
shift evolution of FSRQs have well been developed by Aiello et al. (2012) based on the 2 years 
survey by Fermi-LAI. The all sky (i.e., An integrated) diffuse gamma-ray emission produced by 
FSRQs is determined to be 


= 1.45 x 10 5 GeVcm 


— 2„—1 


( 2 ) 


in the energy range of 0.1-100 GeV (see more details for the derivation in the appendix). 


According to the flat spectrum IceCube detected (eq. 1), the diffuse single-flavor neutrino flux 
in the energy range of 20 TeV-2 PeV is 

r 2PeV 

J u = An / E v $ V } C dE v = 47rln(10 2 )^$ I/ ,i C = 5.5 x lO^GeVcnrV 1 . (3) 

J 20TeV 

If FSRQs can account for the IceCube detected flux, the ratio of neutrino (20 TeV-2 PeV) to 7 -ray 
(0.1-100 GeV) flux for FSRQs is then 


J v /J 1 = 3.79 x 10~ 2 . (4) 

Applying this ratio to individual FSRQs with gamma-ray flux measured we can estimate their 
neutrino flux. 


Analyzing the 4-year data, the IceCube presented results of searching for astrophysical neu¬ 


trinos from candidate sources (lAartsen et al.l 2014d) . The non-detection only gives upper lim¬ 
its to the candidate sources based on the IceCube sensitivity. For an E~ 2 spectrum, the me¬ 
dian sensitivity is ~ 10 - 12 TeV - 1 cm - 2 s -1 in the northern sky for energies^of ITeV-lPeV, and 
l cm - 2 s -1 in the southern sky for energies of 0.1-100PeV (lAartsen et al.l 2014c). 


10- 11 TeV- 1 —- 2 -" 1 


Thus we pay more attention to the candidate objects in the northern sky since they are more strongly 
constrained by IceCube observations. We will compare the predicted neutrino fluxes to the limits 
given by IceCube for individual and stacked FSRQs, thus examine the FSRQ origin of the diffuse 
neutrino emission. 


*We can rule out that in the hadronic model the gamma-ray emission from FSRQs can be produced by secondary 
leptons and photons from photopoin process, given the result of eq. (4) below. See the discussion in the last section. 
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3. Data and results 


3.1. Candidate sources 


IceCube has searched for neutrino emission from 44 candidate sources, among which there are 
11 FSRQs, 5 in the northern sky and 6 in the southern sky ( Aartsen et al.ll2014c ). The non-detection 
by IceCube puts upper limits of neutrino flux to these sources. With the gamma-ray/neutrino flux 
ratio ofEq.S we predicted the neutrino flux from them using gamma-ray flux for individual FSRQs 
from the Fermi -LAT 2-year point source catalog (2 FGlH The results are presented in Table 1, and 
the comparison between predicted flux and the IceCube limit is shown in Fig 1. 

The sensitivity of IceCube for point sources depends strongly on the source declination, as 
shown in Fig. 11 in A artsen et al . ( 2014ch (also shown in Fig 2 here). The typical sensitivity 
in the southern sky, ~ 10 _11 TeV~ 1 cm~ 2 s _1 , is much larger than the predicted neutrino fluxes. 
However, the sensitivities in the northern sky is much higher, and 3 of the 5 northern sources are 
with predicted fluxes above the upper limits. In details, the predicted flux of the source 3C 454.3 
is larger than the upper limit by a factor of 7.3, PKS 1506+106 by 1.4, and 3C 273 by 2.0. 


3.2. Stacking sample 


Since it is usually null result for detection of individual sources, IceCube also carried out 
stacking searches for different classes of objects (A artsen et aL 2014g). For FSRQs, they select 
33 sources based on the measured Fermi -LAT gamma-ray flux, assuming prevalence of photo- 
hadronic neutrino production. The upper limit in the 90% confidence level, for a E 2 flux normal¬ 
ization of flux, is 3.46 x 10” 12 TeV cm _2 s _1 for the total flux of the whole sample in the 

stacking analysis. So the limit for each FSRQ on average obtained from stacking is smaller by 33, 


Efe „ uni = 1.05 x 10 13 TeV cm 


— 2„—1 


(5) 


This is already smaller than all the predicted flux for individual FSRQs in the sample list (see Table 

2 ). 

Since the detection of muon and anti-muon neutrinos by IceCube depends sensitively on the 
declination of the source, we may focus more on the directions of the sky with higher sensitivity, 
especially the northern sky. The FSRQ sample in the stacking search by IceCube includes 19 in 
the northern sky and 14 in the southern sky. We here derive the upper limit for the 19 northern 
ones. 


^http://fermi.gsfc.nasa.gov/ssc/data/access/lat/2yr_catalog 
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Fig. 1.— The predicted fluxes and the upper limits for candidate FSRQs in Table 1 as function 
of sin(<5), with 5 the FSRQ declination. The circles show the predicted single flavor neutrino flux 
of individual FSRQs, and the inverse triangles the upper limits by IceCube of 4-year running. 
Those marked in red are with predicted flux larger than the IceCube limit: 3C 279, 3C 273, PKS 
1502+106, and 3C 454.3 (from left to right). 
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Table 1: Gamma-ray and neutrino fluxes from candidate FSRQs. 


FSRQs 

r.a. 

dec 

f 7 

,pred 

im 

Northern Sky 

4C 38.41 

248.81 

38.17 

1.44 

0.74 

2.71 

3C 454.3 

343.50 

16.15 

12.02 

6.18 

0.85 

PKS 0528+134 

82.71 

13.55 

0.39 

0.20 

0.80 

PKS 1502+106 

226.10 

10.49 

3.82 

1.96 

1.39 

3C 273 

187.28 

2.04 

2.72 

1.40 

0.72 

Southern Sky 

3C 279 

194.04 

-5.79 

2.98 

1.53 

1.51 

QSO 2022-077 

306.42 

-7.61 

1.08 

0.55 

2.07 

PKS 1406-076 

212.24 

-7.87 

0.24 

0.12 

1.66 

QSO 1730-130 

263.28 

-13.13 

0.41 

0.21 

3.46 

PKS 1622-297 

246.53 

-29.81 

0.30 

0.20 

17.2 

PKS 1454-354 

224.36 

-35.67 

1.21 

0.79 

19.6 


Note - The columns from left to right show: the FSRQ names; the right ascension and declination 
of the sources; the gamma-ray flux measured by Fermi -LAT in the 0.1-100GeV energy range (in 
unit of 10 -10 ergcm -2 s -1 ); the predicted u tJ + flux in unit of 10 -12 TeVcm _2 s -1 , assuming E~ 2 
flux normalization and even flavor ratio of 1:1:1; and the upper limit of ;/ /7 + z/ /( flux (in unit of 
10 l2 TeVcm ~ 2 s~ ') in the 90% confidence level from IceCube, assuming E~ 2 flux normalization. 
The source names in bold face mark those FSRQs with predicted flux violating the upper limit. 
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As the upper limits are basically given by the sensitivity, and the sensitivity as function of 
declination (5) has been given by IceCube (Fig 11 in (Aartsen et al. 2014cj) ). we may estimate the 
upper limit of individual sources by the sensitivity at their declinations. The ratio of the upper limit 


for stacking 19 northern FSRQs, to that for the whole-sky 33 ones, /^j- rn = 33 f I>, y/( j im , is 

19 


j?north ^^,sens {Si) 

lim i =1 


7 ptotal 


f: 


33 


( 6 ) 


iWim y" <f> (A A 

Z_ ^ ^V(x,sens\ u jJ 
3 = 1 


where $ i > M!Sens ((5) is the sensitivity at declination 5, and d t (Sj) is the declination of No. i (j) source. 
The average upper limit for the 19 northern-sky FSRQs, = /“°^/19, is then 


33 

fT)north _ .. 


(V) 


Here q can be evaluated by the 4-year sensitivity of IceCube in 90% confidence level (Fig.2), 
r] — 0.19 (see more details in the appendix). Thus the upper limit of the iq, + iq, flux per source 
from the northern sky is 


= 3.45 x 10 _14 TeVcm _2 s _1 , 


( 8 ) 


which is about 3 times more stringent than the limit from the whole sample of FSRQs (eq 5). 


By the method introduced in Section 2, the predicted + iq, flux per source (i.e., the total 
flux divided by the source number) for the whole sample is 




vet'll, pred 


= 1.17 x 10 12 TeVcm 


(9) 


Note, compared with the limit of eq. (5), the prediction is about 11 times larger. As for the 19 
northem-sky FSRQs, the predicted average flux after stacking is 

^$“”pred = °- 93 x 10 _12 TeVcm _2 s _1 . (10) 

Thus the predicted flux is even 27 times larger than the limit for the northern sky FSRQs (eq. 8 ). 
We show in Fig 2 the predicted fluxes of the 33 FSRQs, as well as the predicted stacked flux and 
the upper limit for the northern-source stacking. 


4. Conclusion and Discussion 

We derive the muon neutrino (and anti-muon neutrino) flux of FSRQs under the assumptions 
that there is a connection between the gamma-ray and neutrino fluxes, and that FSRQs can account 
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Fig. 2.— The predicted fluxes for FSRQs in the stacked sample (Table 2; dots) as function of 
sin(<5), with S the declination. The red solid lines mark the predicted single-flavor neutrino fluxes, 
whereas the blue dashed lines are the 4-year IceCube upper limits. The thin {thick) lines are for the 
stacked sample from the whole (northern) sky. Also shown is the IceCube sensitivity of 4 years of 
running time {black dotted line). 
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for the diffuse neutrino flux detected by IceCube. By comparing the derived flux and the limit 
for FSRQs from 4 years of running time of IceCube, we find that (i) a significant fraction of the 
candidate FSRQs, 3 among 5, has derived flux larger than the upper limit, in particular, 3C 454.3 is 
predicted to be exceeding the IceCube upper limit by a factor of 7; and (ii) in the sample selected 
for IceCube stacking search, the predicted flux is more than one order of magnitude larger than 
the upper limit for the whole sample, and even 27 times larger for the northern sample which has 
stringent limit by IceCube. 

The confliction between the prediction and the observation implies that at least one of the 
assumptions is wrong. If the first assumption is valid, the second one will not hold, i.e., FSRQs 
are not the main sources of the diffuse neutrino emission. By the contrast between the predicted 
flux and the upper limit for the stacked sample, FSRQs produce at most 10% of the IceCube- 
detected diffuse neutrino flux. Furthermore, according to the constraint from the northern-sky 
sample FSRQs can only produce a diffuse neutrino flux < 3.7% of the IceCube-detected one. 

As the sensitivity improves with running time (roughly oc 1 /y/i), the IceCube observation in 
coming years will make more stringent constraint on the FSRQ origin of diffuse neutrino emission. 
In particular, the FSRQs in the sky region that IceCube has better sensitivity are favored to be se¬ 
lected to constrain their neutrino flux, and hence make more stringent constraint on FSRQ neutrino 
flux, e.g., those sources in the region of sin 5 ~ —0.1 — 0.9 according to the IceCube sensitivity 
dependence of declination (Fig 2). 

It should be noted that the sample of 33 FSRQs selected by IceCube collaboration in the 
stacking search is selected with the weight of the measured gamma-ray flux by Fermi- LAT. This 
approach is consistent with our assumption that the gamma-ray and neutrino fluxes are well corre¬ 
lated. So we can use this sample for our purpose. Moreover, we note that if for some reason the 
faint FSRQs produce larger neutrino emission than bright ones, then our method and conclusion 
are not valid. However, since the contrast between the prediction and the limit is by orders of 
magnitude, we may have strong confidence to reject the FSRQ origin scenario for diffuse neutrino 
emission. 

An implication from eq (4) is that the FSRQ gamma-ray emission cannot be produced by the 
secondary leptons and photons in the pion production processes. The energy carried by neutrinos 
and secondary leptons plus photons are comparable, thus one expects comparable neutrino and 
gamma-ray fluxes from pion production processes, J 7 ~ J v , in contrast with eq (4). J 7 3> J v (eq 
4) can only appear in unreal cases where pions loose energy rapidly by radiation before decay into 
neutrinos. 
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A. Calculation of diffuse gamma-ray flux from FSRQs 

The all-sky integrated intensity of the diffuse gamma-ray emission contributed by FSRQs can 
be derived as 



(Al) 


This preprint was prepared with the AAS IATpX macros v5.2. 
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where z is the redshift, T is the source photon index, L 7 is the rest frame 0.1-100 GeV 7 -ray 
luminosity, and F 7 is the gamma-ray flux, 


F~ — 


1 14^1^)’ 


(A2) 


with Dl(z ) the luminosity distance of the source. Following Ajello et al. (120120 the distribution 
of FSRQs can be written as 


d 3 N _ d 2 N w dN ^ dV _ w dN w dR 

dL^dzdT ~ dL^dV X FT X FFF ~ ^ ^ X Ft X FFF’ 


(A3) 


where <f>(L 7 ,z) is the LF, dN/dT is the photon index distribution, assumed to be a Gaussian 
distribution with mean // and dispersion cr and independent of 0 , 


diV 

Ft 


(r-M) 2 

e j 


(A4) 


and dV/dz = 47 t D 2 (z)dD c (z)/dz is the comoving volume element per unit redshift, where 
D c (z ) is the comoving distance. A standard concordance cosmology is assumed, with H 0 = 
71kms~ 1 Mpc” 1 and Om = 1 — (1a = 0.27, thus 


D c (z) 



1 

a/ (^m( 1 + A.) 3 + ao 


chi. 


(A5) 


and D l (z) = (1 + z)D c (z). 


For the LF, we adopt the LDDE model suggested in Aje llo et ah (2012), 


$(L 7 ,z) = $(L 7 ) x e(z,L~f) 


where 


4(L 7 ) 


dJV 

dL 7 


A 

ln( 10 )L 7 



(A 6 ) 


(A7) 


e(z,LF) = 


1 + z 


pi 


+ 


1 -|- z c (^LF) J \ 1 A z c (L 7) 


1 + z 


V 2 


-1 


and 


Z c(Lry) = Z* • (L 7 /10 48 )“. 

Here L 1 and L* are in unit of 10 48 erg s _1 . 


(A 8 ) 

(A9) 


The limits of the integration in Eq. IA1I are L,. min = 10 44 erg s 1 , L 7max = 10 52 erg s 1 , 
z ni in = 0, z max = 6 , r min = 1.8 and r max = 3.0. The parameters in the LDDE model have 






















-13- 


been determined by the fitting of the Fermi- LAT data, A = 3.06(±0.23) x 10 9 Mpc 3 erg 1 s, 
7 i = 0.21 ± 0.12, 7 2 = 1.58 ± 0.27, L* = 0.84 ± 0.49, z* = 1.47 ± 0.16, a = 0.21 ± 0.03, 
p l = -7.35 ± 1.74, p 2 = 6.51 ± 1.97, p = 2.44 ± 0.01, and a = 0.18 ± 0.01. We calculate 
the all-sky diffuse gamma-ray flux from FSRQs using these parameter values, and give J 7 = 
1.45 x 10 -5 GeVcm -2 s -1 . 


B. The FSRQs in the stacking search of IceCube 

We show the sample of the stacking search in Table 2, with the sensitivity of IceCube of the 
relevant declinations shown, as well as the predicted // /t + v jL fluxes. Summing up the sensitivity 
for the northern objects only and for the whole sample gives 

19 

n = ^ -= 0.19. (Bl) 

3 =1 



Table 2: The 33 FSRQs in the IceCube stacking search 


FSRQs 

sin (5) 


,sens 

Northern Sky 




2FGL J1849.4+6706 

0.92 

0.38 

2.30 

2FGL J0957.7+5522 

0.82 

0.53 

1.49 

2FGL J0654.5+5043 

0.77 

0.16 

1.41 

2FGL J1312.8+4828 

0.75 

0.43 

1.40 

2FGL J0920.9+4441 

0.70 

0.49 

1.33 

2FGL J1635.2+3810 

0.62 

0.74 

1.23 

2FGL J0043.7+3426 

0.57 

0.11 

1.20 

2FGL J0719.3+3306 

0.55 

0.36 

1.20 

2FGL J1310.6+3222 

0.54 

0.27 

1.18 

2FGL J1522.1+3144 

0.53 

1.10 

1.17 

2FGLJ 1224.9+2122 

0.36 

1.92 

0.90 

2FGL J0714.0+1933 

0.33 

0.24 

0.93 

2FGL J2203.4+1726 

0.30 

0.33 

0.92 

2FGL J2253.9+1609 

0.28 

6.18 

0.89 

2FGL J0725.3+1426 

0.25 

0.47 

0.87 

2FGL J1504.3+1029 

0.18 

1.96 

0.77 

2FGL J1016.0+0513 

0.09 

0.25 

0.73 

2FGL J1229.1+0202 

0.04 

1.40 

0.64 

2FGL J0217.9+0143 

0.03 

0.33 

0.63 

Southern Sky 




2FGLJ1256.1-0547 

-0.10 

1.53 

1.17 

2FGL J2025.6-0736 

-0.13 

0.55 

1.82 

2FGL J0808.2-0750 

-0.14 

0.70 

1.83 

2FGLJ1512.8-0906 

-0.16 

2.62 

2.19 

2FGL J0730.2-1141 

-0.20 

1.17 

3.16 

2FGL J1833.6-2104 

-0.36 

1.05 

5.53 

2FGL J1923.5-2105 

-0.36 

0.42 

5.53 

2FGL J0457.0-2325 

-0.40 

1.13 

6.28 

2FGL J1625.7-2526 

-0.43 

0.42 

6.96 

2FGL J1246.7-2546 

-0.44 

0.46 

7.07 

2FGL J1457.4-3540 

-0.58 

0.62 

11.64 

2FGLJ1802.6-3940 

-0.64 

0.85 

12.20 

2FGLJ1428.0-4206 

-0.67 

0.67 

13.31 

2FGL J2056.2-4715 

-0.73 

0.52 

14.41 


Note - The last two columns are the predicted // /t + D lt flux and the 90% confidence level sensitivity 
(for muon and anti-muon neutrinos) of the 4 years of IceCube running time at the declination of 
the relevant source. Both are in unit of 10 _12 TeVcm _2 s _1 . 


